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Abstract — The Timing Counter of the MEG (Mu to Electron 
Gamma) experiment is designed to deliver trigger information 
and to accurately measure the timing of the e + in searching 
for the decay fi + — > e + j. It is part of a magnetic spectrometer 
with the fi + decay target in the center. It consists of two sectors 
upstream and downstream the target, each one with two layers: 
the inner one made with scintillating fibers read out by APDs 
for trigger and track reconstruction, the outer one consisting 
in scintillating bars read out by PMTs for trigger and time 
measurement. The design criteria, the obtained performances 
and the commissioning of the detector are presented herein. 

Index Terms — Timing Counter, scintillator counter, high reso- 
lution timing detector. 

I. Introduction 

IN the Minimal Standard Model (SM) the Lepton Flavor 
Violating (LFV) decays are proportional to ( j^) 2 , (where 
m u is the neutrino mass and raw the W boson mass) so 
that the predicted LFV Branching Ratios (BR) are negligible 
« 1(T 40 ). 

On the other hand most theories proposed as SM exten- 
sions (e.g. SUSY) predict much larger LFV BRs. In (1 1 the 
BR(/i + -> e + 7) is predicted in the range 10~ 12 - 10~ 14 that 
is accessible by a high precision experiment, like MEG (Mu 
to Electron and Gamma) Q, Q. 

Observation of this decay would be a definitive proof of 
new physics beyond the Standard Model. The aim of MEG is 
to measure this BR down to a few times 10~ 13 improving 
the current limit of 1.2 x 10~ u estimated by the MEGA 
collaboration @). 

II. Event signature, background and 

EXPERIMENTAL APPARATUS 

The basic requirement for the MEG experiment is an intense 
source of like the one currently available at the Paul 
Scherrer Institute in Villigen (CH), capable of delivering up to 
10 8 /.t + /s in DC mode [5|. The /i + s are stopped in a thin target 
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inside a magnetic spectrometer and decay at rest, hence the 
/i + — > e + 7 signature is given by the two-body kinematic, with 
one photon and one positron simultaneously emitted back- 
to-back with the same momentum, 52.83 MeV/c w m„+/2. 
The recognition of the /i + — > e + 7 relies on the precise 
measurement of the kinematic variables that identify the event: 
the e + and the 7 energy and the relative emission angle and 
timing. 

The major challenge is given by the high /i + stopping 
rate required to collect enough statistic in a reasonable time 
interval, resulting in a large background that has to be care- 
fully identified and rejected. There are two major sources of 
background, in which a positron and a photon can fake a 
true /i + — > e + 7 event: the physical or prompt background 
originating from radiative decay ji + — > e + vv^/ with the 
neutrinos carrying small energy, and the accidental background 
arising from coincidences between a high energy e + from 
Michel decay /i + — > e + vv and a high energy photon from 
other sources like radiative decay, positron annihilation in 
flight or bremsstrahlung. In the former background source, 
the time of the e + and the photon at the decay vertex are 
coincident as for the signal, hence the relative timing is not a 
variable resolving signal from this background. 

In order to minimize the background on the experiment, a 
detector with excellent spatial, timing and energy resolutions is 
required, as proposed in [6| (all values quoted are Full Width 
at Half Maximum, FWHM): 

• Photon energy: 4%; 

• Positron energy: 0.7—0.9%; 

• Single particle timing (equal for photon and positron): 
100 ps; 

• Photon-Positron relative timing: 150 ps; 

• Photon-Positron relative angle: 17—21 mrad. 

With these performances and the needed fi stop rate (10 7 — 
10 s fi/s) the rate of accidental background is dominant with 
respect to the prompt one; in this scenario a high timing 
resolution helps in rejecting fake events. Overall background 
rate is evaluated to be ~ 3 x 10~ 14 under these conditions. 

fi + decays occur in a 205 /zm thick polyethylene target 
placed at the center of a quasi solenoidal magnetic field 
provided by the COnstant Bending Radius (COBRA) magnet 
|7|. Positrons are at first tracked by a set of Drift CHambers 
(DCH), and then impinge on the Timing Counter (TC) [8 1 that 
provides position and high resolution time information (9). A 
schematic drawing of the detector is shown in Fig. [T] and in 
Fig§ 
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Fig. 1. Schematic sectional view of the MEG experiment displaying the 
main components. 




Fig. 2. 
event. 



3D view of the MEG experiment with a reconstructed candidate 



The particular configuration of the COBRA field enables the 
e + bending radius, projected onto the plane transverse to the 
z axis Q to be mainly determined by the total momentum with 
only a mild dependency on its polar emission angle (9). This 
principle defines an energy threshold for positrons to impinge 
on the TC thus strongly reducing the detector occupancy and 
the trigger rate from positrons out of the signal region. 

Fig(3] shows that the threshold effect is somehow blurred 
by the spread in the decay points in the target and the sec- 
ondary interactions in the detector material, leading to lower 
energy positrons passing the cut and contributing to detector 
crowding. Nevertheless the positron TC efficiency drops of a 
factor 2 (10) for positrons - 5.0 MeV (~ 10.0 MeV) below 
the signal energy. 

Photon energy, timing and direction are reconstructed in the 
Liquid XEnon Calorimeter (LXEC) 1101 . a vessel, located out- 

1 A cylindrical coordinate system with the z axis along the beam is assumed 



side the magnet, containing ~ 900 liter of liquid Xenon whose 
scintillation light is read by 846 Photomultipliers (PMTs). The 
geometrical acceptance of the LXEC covers approximately the 
solid angle delimited in azimuth by 120° < 4> < 240° and 
| cos0| < 0.35 for fi + decaying at the center of the target. 

Most of the detector signals are digitized by fast ADCs 
based on the Domino Ring Sampler (DRS) chip: the longitudi- 
nal TC and the LXEC sampling frequencies are set at 1.6 GHz 
to preserve the time information. The slower signals from the 
DCH are digitized at 800 MHz flj]. The DRS incorporates 
on a single chip an innovative circular sample-and-hold to 
continuously sample the signals. When a trigger is fired the 
sampled waveforms are digitized and stored. 

The trigger system digitizes signals from all detectors using 
a set of 100 MHz, 10-bit flash ADCs and process them 
with dedicated FPGA (Field Programmable Gate Array): this 
flexible scheme allows to implement many decision algorithms 
to quickly select interesting events. The one dedicated to 
candidate /i + — > e + 7 events is based on fast time and direction 
information provided by the longitudinal TC and LXEC PMTs. 

A set of calibrations have been developed for each detector 
to correct for time offsets, space misalignment and energy 
scale. In the following chapter the TC design and calibrations 
will be described in detail, with a summary of the perfor- 
mances during last year run. 
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Fig. 3. Efficiency of the TC versus the e + energy relative to the signal 
energy (52.83 MeV). 



III. The Timing Counter concept and design 

The TC provides a precise measurement of the positron tim- 
ing and impact point, to be used for the event reconstruction, 
and the direction and timing information of the positron, for 



a fast /i H 



7 trigger. This information must be provided 



over the full geometrical acceptance for signal events with the 
photon pointing to the LXEC inner face. 

The angular range of the positron acceptance cannot be 
defined easily in terms of the LXEC geometry because of the 
spread of the decay vertex on the target, in particular along 
z. The angular range to be covered by the TC is therefore 
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determined by Monte Carlo (MC) simulations taking into 
account the mechanical constraints. 

The MC simulations have also been used to determine the 
TC radial position to ensure the full efficiency for the signal 
simultaneously limiting the fraction of positrons from Michel 
decays hitting the TC that contribute to the accidental trigger 
rate. 

The detector consists of 2 identical sectors, placed sym- 
metrically with respect to the target in order to cover the 
corresponding LXEC geometrical acceptance. Each sector is 
divided in two sub-detectors based on fast optical devices 
arranged in two layers: the longitudinal (outer) one segmented 
along the <fi coordinate and the transverse (inner) one measur- 
ing the z coordinate, satisfying different requirements (Fig.Q. 




Fig. 4. Picture of a TC sector: on top the fibers connected to the APD read 
out equipment, below the bars connected to the PMT. 



A. The Longitudinal (<j>) Detector 

The longitudinal detector is made of an array of 15 scin- 
tillating bars (Bicron BC404) with roughly square shape and 
dimensions 4.0 x 4.0 x 79.6 cm 3 ; each bar is read-out by a 
couple of Fine Mesh PMTs (Hamamatsu R5924) glued at the 
ends. The bars are arranged in a barrel-like shape to fit the 
COBRA magnet profile with a 10.5° gap between adjacent 
bars. The choice of scintillator type, PMT and bar sizes has 
been discussed in detail in |[T2l . [13|. 

A schematic picture of the electronic readout is shown in 
Fig. [5] Each PMT signal is passively split in three channels 
with 80%, 10% and 10% amplitude fractions, respectively. The 
highest fraction is sent to the Double Threshold Discriminator 
(DTD) specifically designed for time analysis purposes. 

DTD is a high bandwidth, low noise discriminator with 
two different tunable thresholds. The first threshold is set as 
low as possible, compatibly with the noise level, in order 
to minimize the Time Walk effect (see SubSect [VHI-B| i. The 
second threshold is set to accept only signals associated 
to physical events, rejecting spurious hits with low energy 



deposit. The value of the high threshold determines the de- 
tector efficiency and is related also to the trigger threshold 
as will be discussed in Sect lVIIFDl When the DTD is fired 
a standard NIM waveform is generated and then digitized by 
DRS boards together with the PMT waveform. The NIM and 
PMT waveforms are analyzed offline to extract the hit time. 
Several algorithms can be applied, the default current one is 
based on template fit on the NIM waveform. This method 
allows to reduce the contribution from the DRS jitter to the 
overall time resolution. 

The first 10% copy of each signal is duplicated and fed 
to the trigger system and DRS: in this stage a second active 
splitter is used for level translation needed to cope with the 
input ranges of Trigger and DRS boards; the second one is 
sent to a charge integrator to monitor the PMT aging. 
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Fig. 5. Schematic picture of the TC electronic read-out. Signals from each 
PMT are passively split and fed into dedicated electronics channels. In order to 
optimize time and pulse reconstruction performances different splitting ratios 
have been used. 

The described architecture provides accurate and reliable in- 
formation from both NIM and PMT waveforms for each event 
and their combination offers adequate information valuable to 
further improve overall performances. 

B. The Transverse (z) Detector 

The transverse detector consists of an array of 256 scin- 
tillating Multi-Clad fibers (Saint Gobain BCF20) coupled to 
Avalanche PhotoDiodes (APD Hamamatsu S8664-55), placed 
orthogonally with respect to the scintillating bars, positioned 
on top of them and covering the same surface. 

The small sectional area of the fibers (5x5 mm 2 plus 2 x 
0.5 mm wrapping thickness totaling 6.0 mm in size along z) 
fits the required space resolution and perfectly matches the 
5 x 5 mm 2 APD sensitive area. In order to comply with the 
mechanical constraints, the fiber ends at the APD side are bent 
in two different types: one set has almost straight termination 
while the other is "S" shaped with small curvature radii (~ 
2 cm): this causes a light loss that is recovered by wrapping 
the fibers with a high reflectance film [ 14 1. This results in a 1.0 
mm spacing between fibers to accommodate for the wrapping, 
leading the pitch to 6 mm along z. 

The APD advantages can be resumed as: small size, in- 
sensitivity to magnetic field and fast response. The main 
disadvantage is the low gain, which can reach at most ~ 10 3 
if biased near the breakdown voltage, and high capacitance, 
which makes them sensitive to noise. 

Each APD is closely coupled to a readout electronics 
performing signal amplification and shaping. A single APD 
current pulse is read as voltage across a load resistor and 
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amplified. Bundles of 8 APDs are mounted on distinct boards 
carrying amplifying stages, power supplies, HV bias and 
ancillary control signal. Due to mechanical constraints, the 
fiber bundle mounted on a single front-end board is interleaved 
with the fiber bundle connected to the adjacent board: the 
resulting assembly, 16 consecutive fibers read-out by two 
adjacent boards, covers 9.6 x 40.0 cm 2 . 
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Fig. 7. Distribution of (z/jj, e7 . — Zf, ar ). 



Fig. 6. Example of a block of 2 trigger waveforms from APDs, related to 
interleaved bunch of fibers, and covering a Az 9.6 cm wide region on 
one (DownStream) sector. Each waveform results from the addition of signals 
from 8 APDs; vertical scale is 100 mV/div and horizontal scale is 0.2 /xs/div. 
A 100 mV offset is added to keep the whole signal and baseline into the 
ADC dynamic range 0-1 V. 

The analog signals from the 8 APDs on each board are 
summed and fed to the trigger digitizer: an example is depicted 
in Fig J6j this is a block of 16 interleaved fibers read out by 2 
consecutive boards, whose signals can be summed up in order 
to increase the Signal-to-Noise ratio and consequently the 
trigger reliability. Noise levels are measured to be in the range 
20-30 mV r.m.s. Using the analog summed output the impact 
coordinate z is obtained at the trigger level with 0(3 cm) 
precision. 

The front end boards provide also a digital signal from 
a discriminator for each APD identifying each single fiber 
hit. The encoded information is stored on dedicated FPGA 
based VME boards. The coordinate of the fiber Zfn, er can be 
matched to the information obtained from the time difference 
between opposite PMTs on the same bar of the longitudinal 
detector zt ar - An example of a preliminary Zfn, er — z^ ar 
distribution is shown in Fig|7] Since the trigger was fired 
by cosmic rays hitting the bars, there are some tails due to 
the inefficiency of the APD detector: a O (10%) is known 
to be due to geometrical effects (incomplete coverage due to 
mechanical constraints and 1 mm spacing between fibers) but 
the larger amount comes from inefficient event reconstruction. 
A more detailed study on this topic is under way. 

IV. The Timing Counter in the trigger 

The trigger system is designed to operate a quick and 
highly efficient selection of candidate events, which is strongly 
needed for an experiment looking for ultra-rare decays in 
a huge beam-related background environment. The severe 
requirement of a maximum latency < 450 ns prevents us from 
using any piece of information associated with the DCH, thus 
leaving the TC as the only detector available for the positron 
at the trigger level. The use of the TC in the trigger has a 
twofold objective: PMT and APD signals are digitized by 
means of 8 VME 6U trigger boards (named "Typel", 4 for 
each), each hosting 16 input channels equipped with 10-bit 
flash-ADCs at 100 MHz frequency (for details on the trigger 



hardware and architecture refer to lfT31 ). The output data are 
recorded, and processed, by a Virtex FPGA to operate an 
on-line reconstruction of the time and impact point of the 
impinging positrons. Dedicated real-time algorithms have been 
developed to implement baseline subtraction and amplitude 
calibration to compensate for any PMT gain and bar light yield 
difference. 

A TC cluster is defined in association with a positron loos- 
ing energy in a bar such that the signal amplitudes from each 
PMTs exceed a pre-defined threshold. In addition, the sum 
of the two amplitudes must be higher than another threshold, 
which ensures an event selection almost independent of the hit 
position. That threshold is equivalent to an energy deposit of 
« 2 MeV, which makes the selection fully efficient for signal 
positrons. If that condition is met by more consecutive bars 
(which is the usual case), these are grouped in a single cluster 
where only the pieces of information (time and position) of 
the first to be crossed are used. 

The on-line estimator for the positron time T e + is obtained 
from a parabolic interpolation of the leading edge of the sum 
of the PMT signals on the hit bar. A ps 2.5 ns resolution 
is obtained on that observable, which allows us to define 
a conservative 20-ns coincidence window with the 7 in the 
LXEC detector. 

The ^-coordinate reconstruction proceeds in two parallel 
streams to provide, in combination with the bar index, a 
redundant stereo read-out of the positron hit position on the 
TC: 

• APD signal discrimination (an example of which is shown 
in Fig. [6] as recorded by a trigger board); 

• the use of a look-up table to compute the logarithm of the 
PMT charge ratio which, according to Eq. [6] is linearly 
correlated to the z-coordinate. 

In the former case, due to the fan-in of APD signals (8 APDs 
to 1 trigger input channel), a resolution a z s» 3 cm is expected. 
In the latter, the resolution has been measured to be slightly 
worse, a z s» 5 cm, but still suited to obtain an efficient 
background rejection. 

The selection of e + — 7 pairs proceeds as follows, in the 
absence of DCH pieces of information, we need to assume 
that each positron is a signal-like one: meaning, it must be 
back-to-back with the 7 (whose direction is extracted from 
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the position of the PMT in LXEC with the maximum output) 
and carrying the signal momentum. Under this assumption, 
the association of a hit position on the TC with a photon 
entering the LXEC with a given direction is univocal and it is 
straightforward to build an association between the index of 
the maximum amplitude PMT in LXEC and the TC bar and 
APD indices based on Monte Carlo simulated signal events, 
as shown in Fig|8] This mapping is summarized in a look-up 
table which is implemented on-line to check the back-to-back 
condition. 
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Fig. 8. Expected hit map (z, ip) for positron signals on the TC obtained from 
a Monte Carlo simulation. Different color markers refer to positrons emitted 
at different <f> angles in the detector frame, each band being associated with 
a PMT row in the LXEC. 



V. Timing measurement 

Considering a positron impinging on a TC bar, the time 
measured by the inner (in) and outer (out) PMTs can be 
written as: 



ti„ = T T C + bin + TWii 



tout — Ttc + bout + TWou 



VeS 



V c f[ 



(1) 



where Ttc is the actual time of impact of the positron on 
the bar, b in out are offsets depending on the readout chain and 
fixed for each couple of PMTs, TWi n ,out are contributions 
from Time Walk effect, v c g is the effective velocity of light 
in the bar and L is the bar length; the z axis points along the 
main axis of the bar and its origin is taken in the middle of 
the bar. 

From Eq. [T]the true impact time is given by: 



Ttc — 



tin tout bi n ~T" b Q 



TW n . 



2v, 



h 



The determination of the factors b and TW and their 
stability through the whole data taking period will be discussed 
hereinafter in more details. 

VI. Position measurement 

The positron impact point is obtained from the time differ- 
ence between the two pulses as 



(tin ~ tout - (bin ~ bout) - (TW m - TW out )) (3) 



The impact point can also be evaluated using the ratio 
between the charges delivered at each PMT lfl6l : 



Q m = E ■ G in ■ exp 
Qout = E ■ G out ■ exp 



'off 



(4) 
(5) 
takes 



where E is the energy released inside the bar, Gi n ,out 
into account several contributions (i.e. the scintillator yield, 
PMT quantum efficiency and gain), A c ff is the effective 
attenuation length of the bar. Taking the ratio we obtain: 

Qin G 



which leads to: 
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(6) 



(7) 



Qin Qin 

Moreover, from Eq. |4] the energy release in the bar can be 
estimated independently from z: 
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Fig. 9. Logarithm of the PMT charge ratio vs the z coordinate measurement 
using Eq.[3]for bar number 1. The linear fit returns A e ff and the ratio of PMT 
gains. 
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Note that the combination of Eqs.|7]and[8]provides a way to 
evaluate the ratio between A G ff / v c g for each TC bar. Assuming 
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Ueff = 14.0cm/ns, the value of A e g is extracted from a linear 
fit as shown in Figj9] The A c ffS extracted in 2010 for all bars 
fall in the range 40 — 90 cm. 

Two methods for the impact point reconstruction are used 
in different stages of the data acquisition chain. The on- 
line algorithm for z reconstruction in trigger, requiring fast 
response and moderate precision, relies on the charge ratio 
method, while the offline analysis, aiming at the best possible 
resolution, exploits the PMTs time difference. The need for 
offline calibrations to guarantee the ultimate performance (both 
for z and time resolutions) is satisfied by using several tools, 
presented in the following. 

VII. Beam Test performances 

All TC single modules (bar-PMT assembly) were tested 
using narrow particle beams to evaluate the timing resolution 
under conditions close to the experimental ones. The tests 
were carried on at the Beam Test Facility (BTF) in the INFN 
Frascati National Laboratory. The BTF can provide bunches of 
e~ or e + with energy up to 750 MeV and multiplicity tunable 
in the range 1 - 10 10 particles/bunch H3, lfl8l 

In our test, we used it as a single electron source in 
order to simulate the bar crossing by one particle at once. 
The single bar rate under experimental conditions has been 
evaluated to be O(100 kHz) by means of MC simulations. 
Since the PMT response is as short as 10 ns, there is a 
negligible probability to have overlapping by multiple hits, 
thus we correctly reproduced the actual working conditions. 
Rate effects have been however studied to make sure the 
detector performance is not affected by particle crowding and 
are reported elsewhere [12|. 

The signals from PMTs were fed into the DTD and 
then used to provide Start-Stop signals to a TAC (Time-to- 
Amplitude Converter) system that together with an ADC board 
allowed to obtain a 12.5 ps/bin resolution, calibrated using 
high-precision delay lines. The timing resolution was evaluated 
from AT = distributions taken in various positions, 

where f,* n and t out are defined in Eq. [T] This method is 
almost equivalent to the one making use of impact time from 
Eq. [2] with the following differences: i) AT is not affected 
by Time Walk effect at the leading order approximation and 
ii) a significant contribution to time spread can be due to 
the finite spot size. The Time Walk treatment under experi- 



mental conditions is discussed in VIII-B while the beam spot 
size contribution was evaluated as follows: the impact point 
position was determined with an external telescope defining 
a 5 mm spot on the bar under test. The effective velocity 
of light in the bar was estimated from AT vs z, to be 
v c g ~ 14.0cm/ns. The resulting contribution from the impact 
point uncertainty to the measured timing resolution is given by 
o-(AT) = 0.5/vT2cm/(u eff (cm/ps)) = 0(10 ps) under the 
approximation of a uniform beam spot (worst case), and thus 
negligible. A similar contribution was evaluated for the elec- 
tronic chain (DTD-TDC-ADC) by splitting the signal from one 



PMT and looking to the width of T 



stop 



distribution. 



the impact angle to mimic the experiment conditions. A rear 
detector was used in order to veto multi-particle bunches that 
could affect the time measurement. 

A resume of the measured resolutions is presented in fig. 10 
and in fig. [TT| from which it is clear that the averaged time 
resolution for each bar is o~t < 40 ps. 




Fig. 10. Bar time resolution, as a function of bar number and impact point 
position. 



90 
80 
70 
60 
50 
40 
30 
20 
10 




- 






























• 

- 

























5 10 
Bar number 



A mechanical support allowed scanning each bar by moving 
the impact point along its main axis (z direction) and varying 



Fig. 11. Bar time resolution averaged along the bar length, versus bar number. 



VIII. Calibrations 

The integration of the TC in the MEG apparatus requires 
a full set of calibrations to determine the detector parameters 
necessary to optimize its performances. Hereinafter only those 
calibrations involving the TC alone are described. Additional 
calibrations, involving other sub-detectors (e.g. space and 
time alignments with DCH and LXEC) are described else- 
where Ifl9l 

A. PMT gain equalization 

The PMT gain equalization procedure exploits cosmic rays 
hitting the TC: due of the uniformity of cosmic ray hit 
distribution along the bars the charge and amplitude spectra 
of the inner and outer PMTs are marginally affected by 
geometrical effects. It is therefore possible to tune the PMT 
gains acting on the bias high voltages and equalize the peaks 
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of the Landau distributions (Fig 12 1 within 15% maximum 



(Fig 13 1. A percent level offline equalization is achieved using 
suitable weights for each PMT in the trigger tables. Due to 
the non-solenoidal configuration of the magnetic field, inner 
and outer PMTs are subject to different field intensities and 
orientations and therefore operate in very different working 
conditions. 
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Fig. 12. Example of charge spectrum acquired with a TC bar, superimposed 
with a Landau function fit. 



As shown in Fig.(fl"5j), this function reproduces outstandingly 
the experimental data. 
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Fig. 14. Graphical representation of the Time Walk effect. 
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Fig. 15. Time delay versus amplitude relation with fit superimposed. 

For each PMT, an event-by-event correction is applied 
based on the pulse amplitude and the fitted coefficients. This 
correction depends on the low level threshold (LLT) value and 
the PMT pulse shape. The optimal LLT value balancing Time 
Walk immunity and noise rejection was found evaluating the 
timing resolution of each bar with dedicated tests performed 



Fig. 13. Distribution of the Landau peak amplitudes versus the bar number. ufjing ^ method defjcribed m [p^A 



B. Time walk corrections 

The Time Walk effect is the dependence of threshold 
crossing delay time on pulse amplitude as showed in FigfT4] 
The analog waveforms sampled at 1.6 GHz are averaged and 
interpolated separately for each PMT over many events to 
obtain a template waveform, that can be used to evaluate the 
correction required in time reconstruction: for each template 
pulse the leading edge region is used to produce a plot of the 
delay time versus pulse amplitude dependence that is fitted 
analytically with the function: 

TW = A + ByG+Clogx (9) 

where x is the ratio between the low threshold value and the 
pulse amplitude of the PMT and A, B, C are the fit parameters. 



Fig. [T6] shows the comparison between the time resolutions 
obtained with two different LLT values (10 mV and 25 mV) 
on the double bar sample. The difference between the time 
resolutions with the two threshold showed is systematically 
favoring the higher value of 25 mV. On the basis of a complete 
scan of the time resolution (evaluated on the double bar 
sample) versus the low threshold in 5 mV steps from 5 to 
35 mV the optimal value of LLT was found to be 25 mV. 

C. Timing Counter offset correction 

Inter bar time offsets due to electronic chains are evaluated 
using again cosmic rays. All PMTs are equipped with Huber- 
Siihner Enviroflex_400 signal cables with equal length (10 m) 
and with low loss characteristics ll20l . This minimizes the 
possibility of time offsets drifting during the data taking period 
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Fig. 16. Timing resolutions for different values of low level threshold for 
double bar events. 



(several months) and, together with a continuous monitoring 
of the relevant variables, ensures a stable operation for our 
detector. For each bar, the time difference distribution between 
inner and outer PMTs is acquired. Due to the cosmic rays 
uniformity this distribution is expected to be fiat and centered 
at zero; the mean values of these distributions give a direct 
measure of the relative offsets between PMTs on each scin- 
tillating bar. The time offsets between different bars and TC 
and LXEC detectors are calculated using events from muon 
Michel decays as well as different tools ll2D . 
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distribution is directly proportional to v c g in the corresponding 
bar. 

The monitoring of the stability of the bar offset is achieved 
throughout the full period of data taking with several tools both 
during the data taking itself (i.e. using the Michel positrons 
to check the hitmap) and with dedicated cosmic runs. Full 
calibration of time offsets and TW coefficients is repeated 
before each year-long campaign (to account for electronics 
and detector maintenance, rescaling of gain factors). 

D. DTD calibration 

The background rejection is a crucial point of the experi- 
ment, together with the maximum efficiency for the relevant 
signals: both are related to the fast on-line event reconstruction 
performed at the trigger level. Trigger position resolution is a 
key parameter for direction match with the photon entering the 
LXEC and is evaluated to be a(Az) ~ 7.3 cm by comparing 
the trigger-level reconstruction of impact point and the one 
obtained with calibrated At. 

On the other hand, from the point of view of subsequent 
analysis, TC trigger and DTD conditional efficiency study is 
motivated by the request for all triggered events to fire a DTD 
output in order to allow the best event reconstruction. The 
relevant parameter is the threshold efficiency cdtd defined, 
for each bar, as the ratio between the number of events with 
NIM signals and the total number of triggered events on 
the bar itself, plotted versus the high level threshold. This 
one is set slightly lower than the trigger one, in order to 
avoid acquiring events without NIM pulse information. Fig. 18 
presents the efficiency cdt d versus threshold level. Averaging 
oover all bars we obtain 6dtd > 99.9% for High Threshold 
HLT < 400 mV. 
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Fig. 17. Cosmic hit maps before and after bar time offset subtractions. 
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Fig. 17 shows the effect of accounting for bar offsets on 
cosmic ray two-dimensional hit map. The shape of the detector 
hit map becomes rectangular when the correct time offsets 
(lower panel) are used and the z reconstruction is reliable. This 
is to be compared to the upper panel, where the uncorrected 
time differences are used for the impact point determination 
resulting in a deformed hit map with respect to the actual 
detector shape. Eq. [3] clearly implies that the width of each 



Fig. 18. DTD efficiency vs high threshold value. The efficiency plateau is 
reached for HT < 400 raV . 



IX. Resolutions 

A. Timing resolution 

The timing resolution is extrapolated from multiple bars hit. 
On a sample of two consecutive bars hit by a positron, the time 
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difference AT = T2 — T\ is studied. For each bar the impact 
time is defined as in Eq. [2] The time difference can be written 

as 



track 



(10) 



(where <7j are the single bar time resolutions and Otrack 
the contribution from the track length spread. Under the 
hypothesis that the two bars have similar resolutions o\ ps er 2 
and neglecting the term due to track propagation, the single bar 
resolution is estimated from the spread of the AT distributions 
divided by y/2. More precisely, we perform a fit using a 
Gaussian plus second order polynomial function to take into 
account tails due to spurious hits (less than a few % of 
the total). In order to keep a conservative estimate useful 
also for long-term monitoring, we apply looser cuts with 
respect to the final analysis: for example, by requesting a track 
reconstructed in the DCH, the polynomial tail disappears and 
the core resolution remains unchanged. This method slightly 
overestimate the timing resolution since the track length spread 
term is not corrected for. Using triple bar events and evaluating 
the quantity: 

T\ + n 



AT = T-2 



(11) 



the effect of the different path length between bars can be 
rejected at the first order, resulting in a better estimate of 



timing resolution, as showed in Fig. 19 where the red and 
black markers represent respectively the resolutions obtained 
in triple and double bar samples. However the triple sample 
has significantly smaller statistics so the double bar sample 
is kept as the reference tool for checking the detector per- 



formances. The resolutions reported in Fig. [T9] are obtained 

this 
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from a different sample than those shown in Fig. 
was because of the insufficient number of triple bar events 
in the latter one. However the two data sets show very similar 
behavior pointing out a stable operation of the detector. 

A certain degradation of the timing performance between 
the Beam Test (at = 40 — 50 ps) and the Physics Run 
configuration (a t = 60 — 80 ps) has been observed. A few 
factors have contributed to this effect: the need for a lower 
PMT gain in order to withstand the high rate and match the 
dynamic range of the DAQ/electronics chain; a slightly higher 
value for the low threshold due to additional noise from the 
surrounding environment; and intrinsic uncertainties in both 
double and triple bar estimates. A deeper study on alternative 
schemes for extracting the time information is under way in 
order to improve the detector performance. 

B. Impact point resolution 

From the triple bar events it is also possible to estimate the 
resolution on the z coordinate of the impact point, substituting 
the reconstructed time with the reconstructed z coordinate: 



Az 



Z2 



Zl + Z 3 



(12) 



where i = 1,2,3 is the bar index. The position resolution 
obtained from the triple bar sample is 1.1cm < <j z < 2.0 cm. 
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Fig. 19. Timing resolution on double (black markers) and triple (red markers) 
bar events, with LLT=25 mV. Due to the different rates for two and three-bar 
events, the sample used here is different than the one in Fig. |16| 
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Fig 20 shows an example of a Az distribution. 



Fig. 20. z resolution on triple bar events with a Gaussian fit superimposed. 



X. Conclusion 

The TC is a cutting edge detector in the field of high 
resolution timing measurements. It demonstrated a timing 
resolution as good as 40 ps in dedicated beam tests and 
60 ps < a(t) < 80 ps under real conditions in the MEG 
experiment together with a precise measurement of the impact 
point with a(z) < 2 cm. The TC fast response satisfies all the 
requirements of a trigger designed to select events on-line in 
a huge background. The TC has been commissioned and is 
currently operating as part of the MEG detector. 
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